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Species distribution models (SDM) are geographical mod¬ 
els of biospatial patterns in association with environmental 
factors (Franklin 1995). The predictive models of species’ 
geographic distributions are important for a variety of appli¬ 
cations in ecology and conservation (Graham et al. 2004). 
Therefore, species’ distributions can be modelled to provide 
suitable information for many rare and poorly known taxa. 
Some snake species are particularly difficult to detect due 
to their low densities, elusiveness, or long periods of inac¬ 
tivity (Seigel 1993). Thus, their distribution ranges may be 
underestimated and less well known than in other reptiles 
(Santos et al. 2006, Bombi et al. 2009). SDMs can be used 
to fill these knowledge gaps by mapping potential distribu¬ 
tion ranges and so identify sites at which searches are more 
promising than at others and should be considered for con¬ 
servation programmes (Peterson et al. 2000). 

The genus Spalerosophis Jan, 1865 (type species Spalero¬ 
sophis microlepis) includes she species, S. arenarius (Bou- 
LENGER, 1890), S. atriceps (Fischer, 1885), S. diadema 
(ScHLEGEL, 1837), S. dolichospilus (Werner, 1923), S.joseph- 
scorteccii Lanza, 1964, and S. microlepis Jan, 1865 (Sindaco 
et al. 2013, Uetz 2015). The genus occurs in arid and sem- 
iarid regions, the Saharo-Sindian region, from North Af¬ 
rica in the west through Arabia, Iran, Pakistan, to central 
India in the east (Baig & Masroor 2008, Sindaco et al. 
2013, Uetz 2015). Four taxa of two species of the genus have 
been recorded from Iran, i.e., Spalerosophis diadema cliffor- 
dii (ScHLEGEL, 1837), Spalerosophis d. diadema (Schlegel, 
1837), Spalerosophis d. schirazianus (Jan, 1863) and Spale¬ 
rosophis microlepis Jan, 1865 (Marx 1959, Latifi 2000, Fi- 


Rouz 2005, Baig & Masroor 2008, Rastegar-Pouyani et 
al. 2008, ScHATTi et al. 2009, Schatti et al. 2010). The lat¬ 
ter species is distinguished from Spalerosophis diadema by 
distinctive morphological characters including 41-45 mid¬ 
body scale rows (Marx 1959, Schatti et al. 2009; Fig. 1). 

However, Spalerosophis microlepis is a rare and poor¬ 
ly known species and its range is uncertain (Marx 1959, 
Gholamifard 2011). It occurs in western and central Iran 
as well as the Zagros Mountains, in Ham, Lorestan, Fars, 
Khuzestan, Hamadan, Markazi, Qom, Kerman, Ghahar 
MahalTva-Bakhtiyari, Kohkiluyeh-va-Boyer Ahmad, and 
Isfahan provinces (Latifi 2000, Rastegar-Pouyani et 
al. 2008, Gholamifard 2011, Moradi et al. 2013, Kaze¬ 
mi et al. 2015). Records of this species from Semnan, west¬ 
ern Yazd and northern Hormozgan need to be confirmed 
(Baig & Masroor 2008, Schatti et al. 2009). Addition¬ 
ally, the species might be present in Iraq, although this 
requires confirmation (Firouz 2005, Baig & Masroor 
2008, Schatti et al. 2009, Gholamifard 2011). Accord¬ 
ing to Latifi (2000), S. microlepis has been reported to oc¬ 
cur in mountainous areas, foothills, fields, grasslands, and 
semi-desert regions. 

The aims of this study are to provide a comprehensive 
distribution map of S. microlepis, to confirm the presence 
of S. microlepis in doubtful localities, and to identify the 
environmental variables associated with the predicted dis¬ 
tribution of S. microlepis using a maximum entropy distri¬ 
bution modelling approach. 

All records of S. microlepis are based on our own field¬ 
work as well as those from the literature (Frynta et al. 1997, 


© 2017 Deutsche Gesellschaft fiir Herpetologie und Terrarienkunde e.V. (DGHT), Mannheim, Germany 
Available online at http://www.salamandra-journal.com 



Correspondence 


Table 1. Percentages of contributions of variables included in the best-fitting distribution model for Spalerosophis microlepis. 


Environmental variables 

Percent contribution 

Permutation importance 

Biol8, precipitation in the coldest quarter 

40 

38.1 

Biol2, annual precipitation 

20 

4.4 

Bio8, mean temperature in the wettest quarter 

15.6 

0.4 

Slope 

10.5 

23.6 

Bio 17, precipitation in the driest quarter 

8.5 

0.8 

Bio2, mean diel temperature range (monthly mean [max.-min.]) 

3.1 

0.4 

Bio 14, precipitation in the driest month 

1.1 

0.3 

Bio5, maximum temperature in the warmest month 

0.6 

16.4 

Bio7, annual temperature range 

0.5 

15.6 

Biol5, precipitation seasonality 

<0.1% 

<0.1% 


Latifi 2000, ScHATTi et al. 2009, SiNDACO et al. 2013). Ad¬ 
ditionally, we included point localities based on museum 
specimens. The records that were used in this study are from 
the following museums: The Natural History Museum, Lon¬ 
don (BMNH); Field Museum of Natural History, Chicago 
(FMNH); Museum d’Histoire naturelle, Geneve (MHNG); 
National Museum of Natural History, Washington D.C. 
(USNM); Museo ed Istituto di Zoologia Sistematica delF 
Universita, Torino (MZUT); Zoological Museum Shahid 
Bahonar University of Kerman (ZMSBUK); Razi Universi¬ 
ty Zoological Museum (RUZM); Zagros Herpetological In¬ 
stitute Museum (ZHIM); and Department of the Environ¬ 
ment of Qom Zoological Museum (DOEQZM). A total of 
33 locality records for S. microlepis were gathered and used 
in the maximum entropy distribution modelling approach 
(Maxent). 20 environmental variables, describing tempera¬ 
ture, precipitation, seasonality, altitude, all with 30-arc-sec¬ 
onds resolution, were obtained from the Worldclim data set 
(http://www.worldclim.org/; Hijmans et al. 2005). In addi¬ 
tion, a slope layer was built from altitude layer information 
in ArcGIS 10 using the spatial analyst toolbox. First, corre¬ 
lations between all 21 environmental variables were meas¬ 
ured with Pearsons correlation coefficient in SPSS 16. The 
variables with a correlation coefficient > 0.75 were excluded 



Figure 1. Male specimen of Spalerosophis microlepis from central 
Iran, Qom. 


from species distribution modelling (Rissler et al. 2006). 
Then, 10 out of 21 environmental variables were chosen and 
used in this study; see Table 1 for more details. 

Maxent is a modeller approach associated only with 
species presence data that enables the construction of well- 
fitted predictive performance and ecological data. It is con¬ 
sidered one of the most efficient approaches for predicting 
species distribution models based on presence data (Elith 
et al. 2006, Phillips et al. 2006, Elith et al. 2011). How¬ 
ever, testing is required to assess the predictive perform¬ 
ance of the model. Therefore, the most usual approach is to 
divide data into ‘training’ and ‘test’ datasets, thus creating 
relatively independent data for model testing (Fielding & 
Bell 1997, Guisan et al. 2006). Gonsequently, Maxent was 
used with default settings when separating records into 
training and test samples (75 and 25%, respectively) with 
ten replicates, which is a technique that has been proven 
to achieve high predictive accuracy (Phillips & Dudi'k 
2008). Gonvergence threshold and maximum number of 
iterations were carried out by default (0.00001 and 500, re¬ 
spectively). We used cross-validation to evaluate the pre¬ 
dictive performance of the model. Jackknife testing was 
used to produce estimates of the average contribution and 
response of each variable to the model. 

Our model was tested with ‘area’ under the receiver- 
operating characteristic curve (AUG) that has been used 
extensively in assaying species’ distribution models, and 
measures the ability of a model to differentiate between 
sites where a species is ‘present’ versus ‘absent’ (Phillips 
et al. 2006, Elith et al. 2006). Models with AUG = 0.5 in¬ 
dicate a performance equivalent to random; AUG > 0.7 in¬ 
dicates useful performance, AUG > 0.8 indicates good per¬ 
formance, and AUG > 0.9 indicates excellent performance 
(Manel et al. 2001). 

The variables that contribute the most include: bioi8 
(40%), bioi2 (20%), bio8 (15.6%), slope (10.5%), bio 17 
(8.5%), bio2 (3.1%), bioiq (1.1%), bio5 (0.6%), bio7 (0.5%), 
and bioi5 (< 0.1%) (Table 1). The AUG value of our model 
was 0.986 ± 0.005. 

Modelling of the potential distribution of S. microlepis 
reveals the most suitable habitat to lie in mountainous re¬ 
gions, including the Zagros highland and northern and 
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southern Afghanistan, which corresponds to the Hindu 
Kush Mountains and northern Syria (Fig. 2). The doubt¬ 
ful records are not congruent with habitat suitable for S. 
microlepis (Fig. 2). The environmental variables with the 
highest gains for S. microlepis are bioi8, bioiy, bio 12, and 
bio8 (Fig. 3) as they are those that will decrease the model’s 
gain the most when they are omitted; this means that these 
variables have a significant amount of information that is 
not represented by the other variables. 

Our results from modelling are highly compatible with 
the known distribution of S. microlepis, with the exception 
of predicted suitability in Afghanistan and northern Syria 
where the species obviously is absent. However, the close¬ 
ly related species S. diadema, which probably has similar 
ecological traits and habitat preferences, occurs there. The 


term niche conservatism is used to describe the tendency 
of species to have similar ecological needs over evolution¬ 
ary time-scales (Peterson et al. 1999, Wiens & Graham 
2005). According to Acevedo et al. (2014), ecological data 
suggests that niche conservatism may be explained by the 
fragmentation in the distribution range of a species’ ances¬ 
tor, which may have been the propellant of the initial stages 
of divergence, without changes of the environmental niche 
of the allopatric populations. On the other hand, predicted 
suitable areas of S. microlepis in Afghanistan and Syria are 
likely not inhabited by the species due to the lack of acces¬ 
sibility in a biogeographical sense. The suitable areas in Za¬ 
gros Mountains are not connected by suitable habitat to the 
highlands in Afghanistan and Syria. Therefore, the species 
could not colonize these areas. 
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Figure 2. Potential distribution of Spalerosophis microlepis resulting from the best-fitting Maxent model. Predicted occurrence from 
low likelihood (white, 0.0) through green, orange to red (1.0). The question marks refer to doubtful records of Spalerosophis microlepis. 



regularized training gain 

Figure 3. Results of Jackknife evaluations of importance of the variables used for our Spalerosophis microlepis Maxent model. 
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The model obtained suggests suitability for occupation 
to be the highest along the Zagros Mountains in western 
Iran, where most records originate. As already mentioned, 
doubtful records such as Semnan, western Yazd, and north¬ 
ern Hormozgan probably do not refer to S. microlepis and 
probably are based on misidentified S. diadema. In addi¬ 
tion, the results of Maxent modelling did not show highly 
suitable habitat for S. microlepis in Iraq, but isolated popu¬ 
lations of S. microlepis probably are located in the moun¬ 
tainous areas of the Kurdistan region, northwestern Iraq, 
which are considered part of the Zagros Mountains and 
known to harbour many species of reptiles and amphi¬ 
bians also present in the Iranian part of the Zagros. A re¬ 
cent study confirms the occurrence of an isolated popu¬ 
lation of the species in northwestern Iraq (Afrasiab & 
Mohamad 2014). Topographically, the Zagros Mountains 
form a barrier between the central plateau and the Meso¬ 
potamian lowlands and a corridor for the southward dis¬ 
tribution of northern faunal elements (Fisher 1968). The 
mountains host the highest number of endemics on the Ira¬ 
nian Plateau and also are considered part of the 20* global 
hotspot region, the so-called Irano-Anatolian biodiversity 
hotspot (Mittermeier et al. 2004, Gholamifard 2011). 
Iran is home to ten endemic species of snakes and seven of 
these occur in the Zagros Mountains: Xerotyphlops wilso- 
ni, Hierophis andreanus, Eirenis rechingeri, Spalerosophis 
microlepis, Telescopus tesselatus. Pseudocerastes urarach- 
noides, and Eirenis (Pediophis) punctatolineatus condoni. 
Consequently, the Zagros Mountains play a major role in 
the separation, isolation, and speciation of the Iranian her- 
petofauna (e.g., Anderson 1968, Rastegar-Pouyani et al. 
2010, Hosseinzadeh et al. 2014). 

Amongst the ten environmental variables that were used 
in this study, the most important factors were precipita¬ 
tion in the coldest quarter (bioi8) and annual precipitation 
(bioi2), as these variables contributed 40 and 20%, respec¬ 
tively. On the other hand, the AUC value of the full model 
was excellent and the standard deviation (SD) of the model 
was also very low, which implied a good performance of 
the model (Manel et al. 2001). Therefore, the results of our 
modelling indicated a trend for S. microlepis of preferably 
selecting relatively humid habitats (mountainous regions: 
Zagros Mountains, Hindu Kush Mountains). However, for 
S. microlepis, precipitation in the coldest quarter (bioi8), 
annual precipitation (bioi2), and mean temperature in the 
wettest quarter (bio8) were the three factors that were more 
significantly associated with its distribution. Commonly, 
environmental variables such as precipitation or temper¬ 
ature are responsible for the distribution patterns exhib¬ 
ited by many reptile species (e.g.. Real et al. 1997, Brito 
et al. 1999, CuiSAN & Hofer 2003, RodrIguez et al. 2005, 
Hosseinzadeh et al. 2014). As has already been shown by 
many authors, reptiles and amphibians are ectotherms and 
absolutely depend on ambient warmth to raise their body 
temperature and then become active, so that they often 
have limited climatic tolerances and are strongly depend¬ 
ent on climatic conditions (Buckley et al. 2010, Luo et al. 
2012, Hosseinzadeh et al. 2014). 


We conclude that precipitation, temperature, and slope 
play the most important roles in predicting the distribu¬ 
tion of S. microlepis as these factors contributed about 85% 
of the environmental factors to the full model. More field¬ 
work is needed throughout Iran and Iraq to shed more 
light on the remaining ambiguities of the distribution of 
S. microlepis. 
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